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Abstract 

A linear transmitter with correlated Gaussian wliite additive and multiplicative noises and 
a periodic signal coupled either additively or multiplicatively is considered. The correla- 
tions have a weak destructive effect in case of the additive signal and a strong constructive 
effect in case of the multiplicative one. Analytical results for the linear transmitter with the 
additive signal agree with those obtained previously by a different approach. We also show 
how to analytically calculate certain expectation values involving exponentials of Gaussian 
processes. 
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1 Introduction 



The fact that noise can play a constructive role in many physical systems is now 
widely recognized. The best-known examples of such phenomena are the stochas- 
tic resonance (SR) [1] and Brownian ratchets [2]. SR is a phenomenon in which 
the response of a dynamical system is optimized by the presence of a specific level 
of noise and has been detected in so many seemingly different systems that it has 
been claimed to be "an inherent property of rate-modulated series of events" [3]. 
However, it has been recently suggested that the functioning of important natural 
devices, e.g., communication and information processing in neural systems or sub- 
threshold signal detection in biological receptors, rely on phase synchronization 
rather than stochastic resonance [4]. The SR is usually measured by the signal-to- 
noise ratio (SNR) for the output signal but new measures are being proposed and 
discussed [4,5,6]. 
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One of the issues that is now debated and that is an origin of many unexpected 
physical phenomena, is the role of correlations between various noises present in 
the system. This point has been particularly stressed in research on Duffing oscilla- 
tor with additive and multiplicative noises [7,8], where correlations are responsible 
for huge changes in the activation rate, on Brownian ratchets [9], where correla- 
tions between the additive and multiplicative noises can lead to current reversal, 
and on a nonlinear rotator [10], where full correlations between two noises amount 
to nullifying one of them, leaving one of the states effectively noise-free and al- 
lowing for the passage of incoming signals of vanishing intensity. Other contexts 
in which correlations between various noise terms appear to be important include 
a quantum dimer under the influence of colored noises [1 1], a coupled neuron net- 
work [12], synchronization of chaotic oscillators [13], detecting the gravitational 
background [14], and medical imaging [15]. The effect of correlations between a 
Gaussian white noise and a Gaussian colored noise has been discussed in Ref. [16]. 

It is sometimes very difficult to find analytical results for nonlinear models and it 
is for that reason that models which admit rigorous solutions are very interesting. 
One of such models is that of an overdamped linear transmitter, which in full form 
reads 



a = (hit) + /2Wei(^))« + fsit) + hitMt) , (1) 

where Ci,2(^) are the noise terms. It should be mentioned that while (1) formally 
leads to a linear equation of motion with time-dependent coefficients, a multiplica- 
tive coupling between the noise and the system means a "hidden" nonlinearity: The 
noise is supposed to represent many unobserved degrees of freedom coupled to 
the transmitting process in a nonlinear manner. Although chemical reaction with 
a fluctuating barrier provided the original motivation for this model, with a being 
the concentration of a reagent, the presence of the additive noise, representing the 
thermal bath acting on the system, and an additive deterministic signal possibly 
lead to realizations that admit negative values of a{t). It, therefore, can no longer 
be regarded as the concentration of a chemical species. Nevertheless, systems of 
this type can still describe processes of a paramount physical interest like the elec- 
tric potential across a molecular membrane [17] or provide a "skeletal" model of 
various enzymatic reactions [18]. This model, in less general forms, has also been 
discussed in Refs. [19]. With one exception (see below), possible implications of 
the presence of correlations between the two noise terms in (1) have not been ad- 
dressed so far. Including such correlations is an obvious extension to the existing 
research on models belonging to the general class (1) and is the subject of the 
present paper. 

There is one important exception to the general statement made in the preceding 
paragraph. Berdichevsky and Gitterman in Ref. [20] considered the linear transmit- 
ter with correlated noises and an additively coupled signal. In order to calculate the 
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correlation function of the process a{t) they obtained a whole hierarchy (a), (^a), 
{^"^a) etc, where is a noise term, and needed a procedure to close the hierarchy. 
They resorted to considering asymmetric, dichotomous, exponentially correlated 
noises and obtained results for Gaussian white noises only as a limiting case. In 
this paper we show how to treat the Gaussian case directly; we think this approach 
is important and interesting, even though the behavior predicted in Ref. [20] is cer- 
tainly more "rich" than that for a pure Gaussian white noise. We also slightly differ 
with the authors of that reference in our analysis of final results. In addition, we 
extend our previous discussion of a model with a multiplicative signal [18] to the 
case of correlated noises. We find that correlations play an important role in this 
system. 

This paper is organized as follows: In Section 2 we introduce the model and discuss 
a couple of its formal features. We give analytical results for the expectation value, 
the correlation function, and the power spectrum of the system with an additive 
noise in Section 3, and analytical and numerical results for the system with a signal 
coupled to the system multiplicatively (parametrically) in Section 4. A short dis- 
cussion in Section 5 follows. Mathematical details are presented in the Appendix; 
although delegated to an Appendix, these details represent an important part of this 
paper. 



2 The model 

Consider a system described by the following equation: 

a = -{ka + K^iit) + f{t))a + \Ut) + git) . (2) 

S(^) represent the multiplicative and additive noises, respectively, acting on 
the system, /(t) and g{t) are deterministic external stimulations; in the present pa- 
per we assume that either one of them vanishes identically. The constants k and 
A allow for an independent manipulation of the noises' strengths. The model (2) 
belongs to the general class of (1) but is simplified to keep only its most impor- 
tant features. The multiplicative noise represents fluctuations of a barrier and the 
additive noise is a thermal bath acting on the system. 

The noises acting on the system need not to be independent. On the contrary, it is 
very likely that both the thermal bath and the process affecting the barrier height 

are correlated, as the molecular mechanisms responsible for these noises are not 
independent. It is the objective of the present paper to examine the effect of this 
correlation on the dynamics of the system. To keep the model as simple as possible, 
we assume that the noises ^i(t), ^2{t) form a two dimensional Gaussian process. 
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point-correlated in time: 



= 5{t - 1') 



1 c 

c 1 



(3) 



where c G [—1,1] is the correlation coefficient. We now perform the Cholesky 
decomposition [21] of the matrix in (3) and use it to construct the noises from two 
independent processes: 
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(4) 



Here ^(t), rjit) are two uncorrelated Gaussian white noises (GWN): (^(t)) = 0, 
ivit)) = 0, mm) = {vitW)) = 5{t-t'), mv{t')) = CandaUhlghercor- 
relations factorize. A more general mechanism of describing correlations between 
the noises has been introduced in Ref. [7], but we will show that the dynamics of 
the system does not significantly depend on the detailed mechanism of introducing 
the correlations. 

We further assume that the external stimulations have a simple periodic form. The 
evolution equation (2) thus becomes 



a^-(ka + K^{t) + Am cos{Qt + 0))a + Xc^{t) + AVl - c^rjit) cos(Oi + 0), 

(5) 

where (f) is an initial phase of the signal. According to what we have said above, 
only one of the amplitudes Aj^, does not vanish. 

The equation (5) can be formally solved. This solution can be used to calculate its 
expectation value and the autocorrelation function. Before we do that, we must first 
solve a couple of problems that will appear during these calculations. 

First, these calculations involve certain expectation values that, to our best knowl- 
edge, have not been calculated so far: One has the form (^^{t)e~-ft ^ 

and the other ^C(^i)C(^2)e"-^*'^^^*'^^^*'^'^*'^, where p(-), ?(•) are certain deterministic 
functions and ^(•) is a GWN. We show how to calculate them in the Appendix. 

Second, the correlation functions should be averaged over the initial phase of the 
signals, as otherwise they will not be stationary [17,22,23]. Therefore, we will use 
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correlation functions of the form 



27r 



{{a{h)a{h))) = ^ / {<ti)a{t2)) dcj) . (6) 



The braces (• • • ) on the right-hand side of Eq. (6) stand for the average over real- 
izations of the noises. 

Finally, calculating the power spectra of the processes described by Eq. (5) is our 
ultimate goal. Let a{t) be a solution to the equation (5). We expect that it becomes 
stationary for t ^ oo. The simplest estimate of the power spectrum is the square 
modulus of the Fourier transform of the process. If we were to calculate the power 
spectrum numerically, we would calculate a particular realization of a{t), reject 
several (perhaps many) initial values until the series becomes stationary, calculate 
the power spectrum of the remaining series and finally average over the realization 
of the noises and the initial phase of the signal. We now follow this programme 
using the analytical, not a numerical, solution. Thus, for the power spectrum we 
obtain 



Piui) — lim — lim 



To+T 

To-T 

To+T To+T 

,iu!{ti-t2) 



T-*oo 2T To- 

To-T To-T 



= lim 7^ lim / dh [ dt2 {{a{To + h)a{To + t2))) cosuj{ti - t2). 

T^oo Zl To— »oo J J 
-T -T 

(7) 

Note that in finite times approximants of the right-hand side of Eq. (7) we must 
always have Tq ^ T ^ and the limit Tq — >^ oo must be taken first. After a simple 
change of variables, we obtain 

•IT 

COSUT X 



Piu) = lim 4^ / dr 



T+2T 

^lirn^ J dt((a(To + it-|r)a(To + it + |r))). (8) 
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If the process a{t) becomes stationary in the asymptotic regime, its autocorrelation 
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function depends only on the difference of its arguments: 



lim ((a ft 



\r)a[t + \r)))^C{T). 



(9) 



t— »oo \ \ 



We thus have 



oo 




2 j COS LOT C{t) dr 





(10) 



If the process a{t) is stationary, we immediately have 



C(r) = lim {{a{t — r)a(t))) = lim {{a {t) a {t + r))) 

t— >oo t— >oo 



(11) 



and we can see that Eq. (10) is in the standard Wiener-Khinchin form. 



3 An additive signal 

We now turn to discussing specific forms of the general system (5). We start with 
a system with an additive periodic signal, described by the following equation: 

d = -(ka + K^(t))a + Xc^(t) + AVl - c^v{t) +Acos{nt + (f). (12) 
The equation (12) has a formal solution 



where Oq is the initial value of the process. For simplicity, we assume in what 
follows that Oo = 0. 

For the expectation value of the process we obtain 





(13) 
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{a{t)) 



A cos{Qt + + 0) Xkc 

\kc Ao.o'^i 



(14) 



where tan = —Vt/{ka — 

As we can see, only the part of the additive noise that is correlated with the mul- 
tiplicative noise contributes to the expectation value (14). There is a maximal ad- 
missible level of the multiplicative noise, /tmax = V^2A5^, beyond which the process 
diverges. In the asymptotic regime, t ^ oo, the system displays stationary oscilla- 
tions. The amplitude of these oscillations grows monotonically as the multiplicative 
noise increases from zero to the maximal level. This is a manifestation of a con- 
structive role of the noise and it is essentially the same as in the case without the 
additive noise [17]. There is one important difference, though: Without the corre- 
lated additive noise, {a{t)) oscillates around zero, and with this noise present, it 
oscillates around a value that depends on the correlation level and strengths of both 
the multiplicative and additive noises. If we now calculate the amplitude of the 
oscillations relative to this flat background. 



oscillations amplitude A 2{ka — \n^) 

background \k\c\ - \K^y + ' 

we can see that the relative amplitude formally diverges for \kc — > (there is no 
threshold and the system can transfer the periodic signal even in the absence of 
the noise) and monotonically decreases to zero as the level of the multiplicative 
noise increases towards its maximal value. In other words, if the strength of the 
multiplicative noise is too large, the periodic stimulation gets drowned in the flat 
background. In this respect correlations between the multiplicative and additive 
noises act destructively on the system. 

If we take an average not only over realizations of the noises but also over the phase 
of the incoming signal, all oscillations vanish: 



We now turn to calculating the correlations for the process a(t). The full expres- 
sions are rather long, and therefore we present the results in the asymptotic {t 
oo) regime only and after averaging on both the realizations of the noises and the 
initial phase: 
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a(t-lT)a{t+lr) 



Xkc 



t— >oo 



+ 



cos fir 



2{{ka 



K^)^ + fi2) 



+ 



+ 



(A«:c)2(A;, - W) 



-(A;a-^«2)r 



(17) 



We can see that the correlation function (17) exists only if the multiplicative noise 
level satisfies k < y/k^ = /tmax/ V^- For the multiplicative noise levels in the range 
Kmax/V^ < K < Kmax, the first moment of the process a{t) converges, but the 
second does not. 

The first term in (17) corresponds to the constant shift introduced by this part of 
the additive noise that is correlated to the multiplicative one; it is equal to a^. The 
second describes the oscillations in the correlation function introduced by the exter- 
nal periodic forcing in (5). The remaining terms describe the diffusive background: 
First of them is present even in the absence of the additive noise. The second con- 
tains contributions from both the correlated and uncorrelated parts of the additive 
noise, and only the correlated part contributes to the remaining one. Note that the 
sign of this term is always positive: the presence of correlations raises the diffusive 
background. The correlations act destructively on the system with a additive signal. 

The correlation function (17), leads to, via Eq. (10), to the following power spec- 
trum: 



-Psignal 

where 



Xkc 



^0 = I — T-r, 1 (19a) 



is the power density introduced at zero frequency by the constant shift. 



is the power density associated with the oscillatory term, and 
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-Pback('^) — 



(19c) 



is the power density associated with the diffusive background. We can now calcu- 
late the signal-to-noise ratio (SNR) as the ratio between the power densities trans- 
ferred by the signal and by the diffusive background at the frequency of the signal: 



SNR(«;^ Z, c) = 



signal 



(20) 



C^K^{ka - Jl 



where Z = (X/A)'^ measures the relative strength of the additive noise and the sig- 
nal. Regardless of the values of other parameters, SNR drops to zero for = ka, 
meaning that the total power density comes from the noise and the signal is totally 
undistinguishable from it. For Z = (20) formally diverges for /t — > (recall that 
since there is no threshold, the system transfers the signal even for a vanishing mul- 
tiplicative noise). SNR remains finite for Z > and since the second term in the 
denominator of (20) is nonnegative, SNR(«;^, Z, c) ^ SNR(«;^, 0, •). If the above 
inequality is strong, the presence of the additive noise acts destructively on the sys- 
tem's capabilities to transfer the signal. In any case, SNR decreases monotonically 
as the multiplicative noise increases: there is no stochastic resonance. 

For weak signals (Z ^ 1) the SNR scales as Z~^, which is readily seen in (20). 
Furthermore, \/Z > 0,k : SNR{k'^,Z,c ^ 0) < SNR(k^Z,0): If the additive 
and multiplicative noises are correlated, the system's capabilities to transfers the 
periodic signal are slightly worse than those of the system with uncorrected noises. 
Numerical simulations show that the system can transfer signals as weak as Z < 
100, but for such a weak signal the difference between the c = and |c| = 1 cases 
is scarcely noticeable. 

Finally, we may ask what is the relative "height" of the 5-peak in the spectrum 
associated with the signal to the (5-peak associated with what manifests itself as 
a constant forcing and results from the presence of the correlations. Not surpris- 
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ingly. 



-Psignal _ ^ 4(A:n 2^^)^ 1 \ 

Po ~ (A«;c)2 ■ {ka - + ^12' ^ ^ 

which once again shows that in the presence of the correlations between the mul- 
tiplicative and additive noises, the signal is drowned by the correlated noise as the 
multiplicative noise level increases. Note that the ratio (21) diverges if \kc — 0; this 
reflects the fact that if any of the noises vanishes, or if the noises are uncorrected, 
there is no constant shift in the output, cf. Eq. (14) above. 



4 A multiplicative signal 



We now consider a model described by the following equation: 



d = -{ka + + A cos(nt + (t)))a + \ci{t) + AVl - r]{t) , (22) 

where all quantities have meanings and properties of the corresponding quantities 
in Eq. (12) above. The only difference is that the signal is now coupled parametri- 
cally to the transmitter. We have discussed this model in Ref. [18], but only for the 
uncorrelated case (c = 0), where we have shown that because of the presence of 
the additive noise, the system responds, in terms of its autocorrelation function, to 
the signal even though the average realization of a{t) goes to zero. 



4.1 Analytical results 



A formal solution to Eq. (22) reads 



a{t)^X J 



e ^"^^ *'^exp 



-K J i{t")dt" -A J cos{nt" + (f))dt" 



X 



{cat') + Vl^v{t'))dt', 



(23) 



where we have assumed that the initial value a(0) = for simplicity. As in the pre- 
vious Section, we will use this formal solution in calculating the expectation value 
of the process a{t) and its autocorrelation function. The calculations run similarly 
to those for the additive noise and make use of the expectation values evaluated in 
the Appendix. For the expectation value we obtain 
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(a(t)) = -iA«:c/e-(^-^''')(*-*') 



X 



exp 



2A 1 /I 
-— sin -n{t - t') cos U + -n{t + t'] 



dt' 



(24) 



If the signal is weak, -C 1, we may expand the expectation value (24), not 

averaged over the initial phase, to obtain 



1 + 



A cos(f2t + + 



yl + tan^ 



tan0 



1] 



(25a) 
(25b) 



We can see that in the limit of weak signals, the transmitter asymptotically displays 
a constant shift and oscillations. The amplitude of these oscillations relative to the 
constant background grows steadily as the multiplicative noise level increases to- 
wards Kmax: The oscillations induced by a parametrically coupled signal are not 
drowned in the noisy background. The presence of these oscillations is clearly a 
constructive effect of the correlations as in the uncorrelated case, c = 0, {a{t)) = 0. 

If we further average Eq. (24) over the phase of the incoming signal, we obtain 



im)) 



1 

— Xkc 
2 





27r 



j (i^exp — — sin -VL{t - t') cos U) + -VL{t + t')\ 



jit: 



27r 



2tv 



>exp 



2^ • Ir. / 
— — — sm -\lt cos ( 

Vt 2 



,(26) 



since the integration over runs over the entire period of the integrand and therefore 
the value of the inner integral cannot depend on \Vt{t + 1'). An integral representa- 
tion of the modified Bessel function 



hiz) 



E -N^ -TT h'^'P cos 0) d<P (27) 
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can immediately be recognized [24]: 

Mm = —cXk j e-('=-5-^)*7o sin ^Qt'^ dt' . (28) 



Furthermore, as 



t 

.1,^2-, 



1 /^f^)e-(^-^«^)*^0, (29) 



we may conclude that 

oo 

Mm ^ -l^^c J e-(^''-^''')*7o sin ^nt'^ dt' (30) 



provided that k < k^^x = V^ka', otherwise the process is divergent, which is 
readily seen directly from Eq. (24). Observe that ii A ^ ka — , the expectation 
value (30) formally diverges in the limit Q ^ 0, even though realizations for some 
particular values of the initial phase (p may be convergent. 

For the correlation function we obtain 



(t - \t) a{t + \t))) = Ji + J2 + Ja , (31) 
where 

t 

J, = A2e-('=-^-^)- I ds e'''('^~^'>Io (^^Zi) , (32a) 





(32b) 
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Fig. 1. Numerical power spectra of the process (22) for various input signal frequencies 
and correlations between the noises. The input signal frequency equals ^ = tt/8 (panels (i) 
through (iii)) and $7 = 27r (panels (iv) through (vi)). The multiplicative and additive noises 
are uncorrelated (c = 0) on panels (i), (iv), partially correlated (c = 1/2) on panels (ii), (v), 
and fully correlated (c = 1) on panels (iii), (vi). Other parameters, common for all panels, 
are ka = 1/2, k = A = 1/4, and A = 1. 

and Zi are certain functions involving trigonometric functions of r and the inte- 
gration variables, but not of t. We can see that this correlation function converges 

for K < \fk'a = /^max/ V^- A reasoning similar to that used in (29) shows that the 
correlation function becomes stationary in the limit t — > (X). In the limit 1] ^ 0, the 
correlation function (31) formally diverges \i A > 2{ka — ^2)- Note also that the 
square of the strength of the additive noise A multiplies all terms in (32). The situa- 
tion is slightly paradoxical: without the additive noise, almost all realizations of the 
process a{t) eventually die out, yet as in the output process the power attributed to 
"signal" and to "noise" scale by A^, the SNR does not depend on the additive noise 
level. 



4.2 Numerical results 



Because of the analytical structure of Eq. (32), any further analysis of this corre- 
lation function should be performed numerically. Instead of numerically calculat- 



13 



ing the Fourier transform of Eq. (31), we have directly integrated Eq. (22) by the 
Heun method [25] with a time step h = 1/256, calculated the power spectrum 
for each realization and averaged over 512 realizations of the noises and the ini- 
tial phase. We have used Marsaglia algorithm [26] to generate the GWNs and the 
famous Mersenne Twister [27] as the underlying uniform generator. Selected re- 
sults are presented on Fig. 1. As we can see, if the external frequency is small and 
the two noises are uncorrelated, the process a{t) does not display any oscillations, 
cf. panel (i). If the noises become correlated, distinct sharp peaks corresponding to 
the driving frequency and its higher harmonics appear, cf. panels (ii) and (iii). For 
a large driving frequency, in the uncorrelated case there is a broad peak centered 
at the driving frequency, cf. panel (iv). If the noises get correlated, a sharp line 
appears, superimposed on the broad peak, cf. panels (v) and (vi). 

If there are no correlations, oscillations introduced by a parametrically coupled sig- 
nal are, by their very nature, damped with an effective damping constant ka — \k'^. 
If the driving frequency is very small, with Q,~^ ^ (/;;„— |«;^)~^, the transmit- 
ter is dominated by the damping term and the resulting power spectrum is non- 
oscillatory (Fig. 1, panel (i) above). For larger frequencies a broad peak, charac- 
teristic for damped oscillations, appears in the power spectrum (Fig. 1, panel (iv)). 
The presence of correlations between the additive and multiplicative noises changes 
the situation. Sharp lines, characteristic for undamped oscillations, appear in the 
power spectrum for both small and large driving frequencies. We can see that the 
correlations have a strong constructive effect on the transmitted signal in case of 
a multiplicative coupling. 

The presence of higher harmonics of the driving frequency on panels (ii), (iii) re- 
flects a nonlinear nature of the multiplicative coupling between the transmitter and 
the external signal. 



4.3 The origin of the constructive role of correlations 



To understand, at least qualitatively, why correlations between the additive and mul- 
tiplicative noises reinforce the parametrically coupled signal, observe that Eq. (22) 
is formally equivalent to 



a — — 




ka Acos{nt + 0) + AVl - c2 r]{t) . 



K K 



(33) 
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If we now introduce 



a — a — cX/k , (34) 

we obtain 



a = - [ka + K^{t) + Acos{Qt + (j)))a-—ka-—Acos{Qt+(j)) + \Vl - c^r]{t) . 

K K 

(35) 

The correlated part of the noise has been formally removed from Eq. (35). We 
can now see that the presence of correlations effectively amounts to (i) lowering 
of the amplitude of the additive noise, (ii) introducing a constant additive forcing, 
and (iii) additively coupling to the transmitter a periodic signal with exactly the 
same phase and frequency as the multiplicative one. The amplitude of the additional 
signal is proportional to the product of the correlation coefficient and the strength 
of the additive noise. As we already know from Section 3, an additive periodic 
signal introduces undamped oscillations into the correlation function which, in turn, 
introduce very sharp lines into the power spectrum. Such lines are clearly visible 
on Fig. 1. 

Note that a transformation analogous to (34), when applied to the system with 
a purely additive signal (12), does not affect the signal. This observation, perhaps, 
clarifies why the presence of correlations in the additive case does not lead to any 
constructive effects. 



5 Discussion 



In this paper we have discussed the effects of correlations between additive and 
multiplicative noises in a model of a linear transmitter (5). We have shown that 
the effects of correlations is different depending on whether the signal is coupled 
additively or multiplicatively. In the former case, the correlations act destructively: 
The amplitude of the outgoing signal increases with the level of the multiplicative 
noise, but the constant shift introduced by the correlations increases even faster, 
the amplitude of the oscillations relative to the shift decreases, and the signal is 
eventually drowned by the noise. Note that the correlations act here against the 
noise: The multiplicative noise enhances the signal, but the correlations diminish 
it. Similar effects can be seen in the power spectrum and the SNR of the outgoing 
signal. The fact that correlations reduce the destructive effects of noise has been 
reported previously; papers quoted in the Introduction provide examples of such 
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systems. What we have here is an example of correlations reducing the constructive 
effect of noise. 

These formal results, and the correlation function (17) in particular, agree with 
those reported in Ref. [20], where the Gaussian white noise has been obtained as 
a limiting case of a dichotomous, colored noise. This should come as no surprise: 
There is a generally held belief (see e.g. [28]) that in a well defined limit the two, 
in principle very different, kinds of noise have identical effect on the dynamics of 
a system. A comparison of our present results with those of Ref. [20] provides yet 
another confirmation of this statement. 

The situation is entirely different for the case of a multiplicative signal. We have 
shown that for a weak incoming signal not averaged over the phase, the outgoing 
signal displays oscillations that do not vanish as the level of the multiplicative noise 
increases. Numerical evidence suggests that without the correlations, the power 
spectrum averaged over both realizations of the noise and the initial phase does not 
display any signs of oscillations for low frequencies and only damped oscillations 
for larger frequencies. This is easy to understand: A multiplicative signal is, by its 
very nature, damped, and if the wavelength of the signal is large, the signal does not 
recover before if gets completely washed away by the damping. The presence of 
correlations changes this picture altogether: The power spectrum displays distinct, 
sharp lines for all frequencies of the input signal. We have showed analytically that 
correlations between the multiplicative and additive noises are effectively equiva- 
lent to adding an extra signal, with exactly the same phase and frequency as the 
multiplicative one, to the system. The amplitude of this extra signal is proportional 
to the correlation coefficient and it is this signal that produces the sharp lines {5- 
peaks) in the power spectrum. 

We have assumed in Section 2 that the two noises present in the system together 
form a two-dimensional Gaussian process. On a formal level, this property is re- 
flected in the relation between the amplitudes of the correlated and uncorrelated 
parts of the additive noise. If we relax this assumption, instead of (5) we obtain 



a= -[ka + n^it) + cos{flt + 4>))a + A^C(t) + \rj{t) + cos{flt + 4>), (36) 

where no special dependence between and is assumed. The results reported in 
this paper can now be recovered if is replaced by A|/(A| -|- A^) and Z = (X/A^y 
is replaced by ( A| -I- A^)/A^. Thus, results reported do not depend on the assumption 
that the joint distribution of the noises is Gaussian. They do depend, however, on 
the assumption that each individual noise is Gaussian and point-correlated in time. 

1 acknowledge a very helpful discussion with Prof. Peter Hanggi that helped me to 
clarify certain ambiguities present in a draft version of this paper. 



16 



A Appendix: The expectation values 



While calculating the expectation value of a{t) and the correlation function of this 
process, one encounters several integrals that involve the expectation value of the 
exponential of an integral of the multiplicative noise ^{t). Specifically, 



exp 



1 \ fl 
j^{mt')dt' )=exp -jW)fdt' 



(A.1) 



where T is a certain time and i^{t) is a certain function, not necessarily continuous 
(cf. [17,29]). Similarly, 



77(^1) exp 



r]i'ti)r]{'t2) exp 



J <p{mt')dt' 


T 



=0, 



5{ti—t2) exp 



\ jW)?dt' 



(A.2a) 



(A.2b) 



because r]{t), ^{t) are independent GWNs. The other two expectation values are 
more challenging. We start with 



Ci = (e(ii)exp 



T 



^{t')m dt' 



(A.3) 



We expand the exponential function in a Taylor series, interchange the orders of 
summations, and obtain 



T T 

°° 1 r r 

Ci = E -r / dt,, ... dt,Mt,^) ■ ■ ■ ^{t,J {^t,,) . . . at^Mti)) ■ (A.4) 

n=o ^- i 



The expectation value on the right-hand side of (A.4) involves n+1 terms. There- 
fore, only the terms with n — 2m+l will give a nonzero contribution to the sum. 

The expectation value thus involves 2(??i+l) noise terms. It factors out to a product 
of m+1 two-point correlations and there are (2m+l)!! ways to choose the pairs. In 
every possible choice, one of the pairs involves ti and a certain t^, and the remain- 
ing m pairs involve two t^'s. Thus 
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m=0 



(2m + 1)! 



T 



\0 

T 



m=0 \ n J 



- 1 [ip{t,)f dt. 



(A.5) 



L 



The last remaining expectation value involves two noise terms multiplying the ex- 
ponential: 




where, as above, we have expanded the exponential into the power series. Only the 
terms with n — 2m will give a nonzero contribution to the sum. The expectation 

value of the product of noise terms factors out to a product of two-point correla- 
tions. Again, there are (2m+l)!! ways to choose the pairs. In (2m— 1)!! cases the 
noises at ti, t2 will be coupled with each other, and in the remaining cases they will 
be coupled with different t^s. Therefore 
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m=0 (2^)' 

T T 







^ g (2m+l)!!- (2m-l)!! ^ 



m=l 



(2m)! 

T \ / 



T T 



m— 1 







00 -j^ 



m—l 



\ 2'"-i(m- 1)! 



\0 



= [S{ti - 12) + 99(^1)99(^2)] exp 



L 



(A.7) 



This approach can be easily generaUzed to colored Gaussian noises. Suppose E,{t) 
is a Gaussian noise with zero mean and correlations of the form 



mm) = h{\t-t'\), 

where h{u)du — 1. Then, instead of (A.7), we obtain 



(A.8) 



exp 



mm dt 



h{\h-t2\)+ j ^it^)hi\t^-h\)dt^j \J ^iU)h{\U - t2\) dU 

\^ T T 

-^jdtnj dU(p{ti^)(p{U)h{\ti^-U\) 



X 



exp 







(A.9) 



and similar expression for other expectation values. 
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As we have said, the function </? needs not to be continuous. In many practical 
situations, for example in evaluating expressions like 



exp 



«2 



(A. 10) 



with ^ tl < t2 ^ T, (p{t) would be the characteristic function of an interval: 






t < tl 


1 

2 


t^tl 


1 


tl <t <t2 


1 

2 


t^t2 





t > t2 



(A.11) 



The rationale for putting \ at the endpoints is that 



1 
2 



t2 OO 

J 6{t-ti)dt= J (p{t)6{t - tl) dt = (p{ti) . (A.12) 

tl —GO 



The first equality stems from the fact that in fl^ 5{t — ti) dt the integrand contains 
only "a half of the peak". By the same token, for example 



^{ti) exp 



J. 

-K J C(t') dt' 
tl 



-lKexp(^lK\T -tl)) , (A.13) 



where T ^ti and Eq. (A.5) has been used. 
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